Linear optics substituting scheme for multi-mode operations 
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We propose a scheme allowing a conditional implementation of suitably truncated general single- 
or multi-mode operators acting on states of traveling optical signal modes. The scheme solely relies 
on single-photon and coherent states and applies beam splitters and zero- and single-photon detec- 
tions. The signal flow of the setup resembles that of a multi-mode quantum teleportation scheme 
thus allowing the individual signal modes to be spatially separated from each other. Some exam- 
ples such as the realization of cross-Kerr nonlinearities, multi-mode mirrors, and the preparation of 
multi-photon entangled states are considered. 
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I. INTRODUCTION 

A general problem in quantum optics is the imple- 
mentation of a denned single- or multi-mode operator by 
some physical scheme since many efforts amount to the 
realization of a desired quantum state transformation of 
an unknown or control of a known signal quantum state 

Here, we limit attention to a physical system consisting 
of a number of traveling light pulses described by non- 
monochromatic signal modes. We assume that the spa- 
tial pulse length is large compared to the wavelength of 
the radiation, so that we deal with quasi-monochromatic 
pulses corresponding to quasi-orthogonal modes. 

A difficulty is that at present, only a limited set of 
basic operations can be implemented directly. Compos- 
ite operations therefore have to be constructed from el- 
ementary ones which are simple enough to allow their 
direct realization. Examples of such basic operations 
are the preparation of coherent (i.e., Glauber) states 
by single-mode lasers, or single-photon states by para- 
metric down converters, furthermore parametric interac- 
tions as realized by beam splitters, three-wave mixers, 
and Kerr-nonlinearities, and the discrimination between 
presence and absence of photons by means of binary 0-1- 
photodetectors such as Avalanche photodiodes. A simple 
combination of these techniques allows further manipula- 
tions such as parametric amplification, coherent displace- 
ment, or the preparation and detection of photon number 
(i.e., Fock) states 

To give an example, the signal pulse may be sent 
through a medium applied in the parametric approxi- 
mation, i.e., the medium realizes a weak coupling be- 
tween the signal modes and a number of auxiliary modes 
prepared in strong coherent states. Since within clas- 
sical optics, nth-order interactions are described as a 
nth-order deviation from linearity of the polarization in- 
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duced in a medium by an electric field, the interaction 
strength is expected to decline rapidly with increasing 
order. Strong fields are therefore required for their ob- 
servation. In contrast, setups discussed within quan- 
tum optics and quantum information processing often 
operate with superpositions of low-excited Fock states 
while the coherent amplitudes in the auxiliary modes 
cannot be increased unlimited to achieve a desired in- 
teraction strength of the reduced signal operation. This 
greatly limits the order of the nonlinearity applicable 
and with it the variety of unitary transformations that 
can be realized by a given medium. Desirable are there- 
fore substituting schemes for such nonlinear interactions 
|l8l Il9| . In particular, one may apply the nonlinear- 
ity hidden in the quantum measurement process and to 
use merely passive optical elements such as beam split- 
ters [2^ and photodetectors while relying on simple aux- 
iliary preparations such as coherent and single-photon 
states 21] . A local measurement performed on a spa- 
tially extended quantum system affects the reduced state 
at the other locations. By repeated measurements and 
postselection of desired detection events, this back-action 
caused by a measurement can then be used to perform 
well-defined manipulations, including state preparations 
lll^BEIJiHjpgjplli, quantum lo gic operations 

state purifications 36] , quan- 
tum error corrections [33 , and state detections |38l l39| . 
The simplest basic operations one may think of are the 
application of the mode operators a and a\ i.e., pho- 
ton subtraction and addition jljj, |4l|. With regard to 
compositions, it is necessary to limit operation to a suit- 
ably chosen finite-dimensional subspace of the system's 
Hilbert space, cf., e.g., [53, since only a finite number of 
parameters can be controlled. 

The aim of this paper is to investigate a theoretical pos- 
sibility of implementing a desired single- or multi-mode 
operation Y on the quantum state g of a traveling opti- 
cal signal. The scheme solely relies on beam splitters as 
well as zero- and single-photon detections. In this way, 
nonlinear and active optical elements can be avoided. It 
requires the preparation of coherent states and single- 
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photon states, however. The idea is to start with a single 
photon, which is then manipulated to construct an en- 
tangled k-mode state. After that, the latter is shared by 
local setups performing the transformation in the signal 
modes. The measurement-assisted and hence conditional 
transformation leads to an output state 



g' = -YgY\ 
P 

where the normalization constant 

p = Tr(YgY^) = (Y^Y) 



(1) 



(2) 



is the probability of the respective measurement result, 
i.e. the 'success probability' . In particular, if the oper- 
ator Y is proportional to a unitary one, Y — all, where 
U> = U~ l , then the success probability is independent of 
the input state, p = \a\ 2 . 

Within the scope of explaining the principle, we limit 
attention to idealized optical devices, as the effect of im- 
perfections depends on the desired transformation and 
the signal state itself. In a given practical setup, the 
imperfections and the resulting coupling of the system 
to its environment must be considered since loss of only 
one photon may change the phase in a superposition of 
two states and yield wrong results. Apart from this, the 
mode matching becomes an important issue especially in 
the case of composite devices such as optical multiports. 
A generalization of Eq. (JJ) allowing the inclusion of loss 
would be the transformation g' = p^ 1 KiYigY^ with 
given coefficients m > 0. 

The article is organized as follows. SectionlTfldescribcs 
the local devices which mix the signal pulses and an en- 
tangled state in order to perform the desired transforma- 
tion of the signal, whereas section ITTT1 is dedicated to the 
preparation of the entangled state itself. The operation 
of the complete setup is considered in section llVl and sec- 
tion[V]cxplains how general operators can be approached 
with it. To give some examples, a few operators of special 
interest are considered in section IVI1 such as functions of 
the single-mode photon number operator in section lYl Al 
two- mode cross-Kerr interactions in section IVl Bl U(A)- 
transformations in section IYl CI and the preparation of 
multi-photon entangled states in section lVIDI Finally, a 
summary and some concluding remarks are given in sec- 
tion lVIII An appendix is added to outline three possibili- 
ties of implementing single-photon cloning as required to 
prepare the entangled state and some ordering relations 
for the photon number operator. 



II. 



LOCAL OPERATIONS 



The scheme consists of k local devices, each implement- 
ing a conditional single-mode photon subtraction or ad- 
dition, depending on the case in which it is applied. The 
setup of such a local device is shown in Fig.^and consists 
of beam splitters Uoi(Ti, R\) and Uzifi^Rz) as well as 



photodetectors Di and D2. We describe the mixing of 
some mode j with some other mode k by a beam split- 
ter in the usual way by a unitary operator Ujj, (T, R, P) 
defined via the relation 



U 



(3) 



Jj Ujk P {-R* T*) [& kJ 

by its complex transmittance T, reflectance R, and phase 
P obeying |T| 2 + |i?| 2 = \P\ 2 = 1. If not stated otherwise, 
we will assume that P = 1 within the work. 

The setup shown in Fig. ^ realizes a transformation of 
the signal state g in mode according to Eq. Q , where 
the single-mode operator Y =Yq(s) depends on the pho- 
ton number s — 0, 1 of the Fock state \s) in which one of 
the remaining input modes are prepared. We distinguish 
between two different cases of operation. First consider 
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Case (1) 



FIG. 1: Setup consisting of beam splitters U and photode- 
tectors D and implementing conditional subtraction [case (0)] 
or addition [case (1)] of a photon on a signal state q, each 
controlled by the photon number s = 0, 1 of the Fock state \s), 
respectively. Note that in case (1) the beam splitter U21 is 
redundant. 

case (0), in which input mode 1 is prepared in the vac- 
uum state |0) and input mode 2 in \s). If Di and D2 
detect and 1 photons, respectively, the operator Y (s) 
becomes 



i?2 



2<l|E/ 2 it/oi|0)i|s) 2 
TV " 

i?2, 



T^(-Rtao) 1 



(4) 



A a 3 



is the photon number operator. In case 
(1), input mode 2 is prepared in the vacuum state |0) 
and input mode 1 in \s). If Di and D2 detect photons 
each, the operator Yq(s) reads 



|£f 2 itfoi|s>i|0> 2 



1? 



77" ( ^4 



(5) 



We see that (apart from the operator T™° which is always 
present) in case (0) we realize a photon subtraction and 
in case (1) a photon addition, controlled by the photon 
number s = 0, 1 of the control input state |s). 
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III. PREPARATION OF THE ENTANGLED 
STATE 

The combined operation of k local setups shown in 
Fig. requires an entangled A:- mode state 



I*) 



• • ■ \s k )k + z\l - si) i ■ ■ ■ |1 - s k ) k 



(6) 



shared by all the stations. Here, z is a given complex 
number, and the \s) are photon number states with s = 
0,1. In this section, we discuss the preparation of the 
state Eq. ©. We start with preparing a mode in a 
single-photon state |1) and mix it with some other mode 
j prepared in the vacuum state |0) using a balanced beam 
splitter U j(T = —R= l/\/2), which leaves the photon in 
a state 



\*) 0j = U 0j \l) \0), 



10)011), + |l)o|0) 3 
V2 



By a repeated application of single-photon doners 

\ s )k\s)j j{s\, 



h 

s=0 



(7) 



(8) 



which duplicate photon number states |s) with s = 0, 1 
according to Qkj\s)j = \s)k\s)j, the state Eq. Q can now 
be enlarged to a (k + l)-mode state 



|*)o k = Qkj- 

|0)o„ 



• Qj+i,jQj-i,o ■ ■ ■ Qio\^)oj 
„ J --i|l)i,...,fe + |l)o,...j-i|0)^..., fe 



V2 



, (9) 



where 0<j<k and k > 0. Note that here we have used the 
notation \s) m ..... n — \s) m \s) rn+ i ■ ■ ■ \s) n . Some possibili- 
ties allowing the implementation of Eq. JSJ) are outlined in 
appendix^ cf . also 42] . To manipulate the state Eq. © 
further, we mix mode with an auxiliary mode fe+1 pre- 

pared in a coherent state \a) =e ~ X)^lo a "( n ')~ 5 l n ) 
using a beam splitter ?7o,fc+i(T', R). If thereafter 1 and 
photons are detected in mode and k + 1, respectively, 
the state Eq. © is reduced to the fc-mode state Eq. ©, 



o<l|fc+i<0|C/o,fc+i|*)o,...,fc|Q!)fc+i 



'v^l*> 



.,fc- 



(10) 



Here, ip — argT is an unrelevant phase and the photon 
numbers si, . . . ,s k can be chosen arbitrarily to be either 
or 1 by permutating the mode indices suitably, except 
that one of the s; is always 0. Note that the latter does 
not constitute a limitation since we have 

1*0,1, • • • ,s k ;z)) = e™|*(l- Sl) • • • , l-sk^z- 1 )). 

The complex parameter z = T 1 Ra can be controlled by 
T or a. It is however convenient to maximize the success 
probability 



(i + H 2 )HV 

2(M 2 + |z| 2 



(12) 



by adjusting \T such that 



4 



2 



(13) 



is valid, in which case it is ensured that (2e) 1 < p < 2 1 
holds for all z. 



IV. 



OVERALL OPERATION 



Let us now assume that a state l*)!,...^ given by 
Eq. © is prepared. Each mode 1 = 1, . . . , k is fed into the 
control input of a device Fig. ^ (the input port fed with 
the variable photon number state | s) in Fig. ^) , assuming 
case (s;), respectively. A schematic view of an example 
of the overall setup is depicted in Fig. Denoting the 
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FIG. 2: Example of the overall setup consisting of k de- 
vices Fig. Q applied in the signal modes ji, . . . , jk in cases 
(si), . . . , (sfc), respectively. These modes do not have to be 
distinct from another. The respective control input ports of 
the devices which are fed with states \s) in Fig. Q here share 
an entangled state |\t)i,...,fc given by Eq. Ir- 
respective signal mode (mode in Fig.^) by ji, we obtain 
a transformation in modes ji , . . . , jk of the signal state g 
according to Eq. (JIJ, where 



Y{A,a,/3) 



Y 



+^)(i- s ,)...y.^)(i- Sl ) 

af(A + p), 



(14) 
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with 



A. Existence of an eigenspace of A 



T 
A 



1 



Ok "ii ' 



E?=i«. 



fc-Ef=i«i 



x2? 



(15a) 
(15b) 

(15c) 



(15d) 



and a^ s ' stands for a if s = and for if s = 1. We see 
that for given beam splitter parameters 71- and Rk , /? can 
be chosen as desired by varying z. Note that the signal 
modes j\, . . . , jk don't have to be distinct, cf. the ex- 
amples in section lYl Al and section lYl Bl Note also that 
those modes which are distinct may be spatially sepa- 
rated from each other. The signal flow of the scheme 
Fig-Elthcn resembles that of a multi-mode quantum tele- 
portation setup. 

We now consider N successive applications of Eq. (jTlfl . 
If z is varied from step to step, the resulting overall op- 
erator becomes the product of the individual operators 
Eq. FP^i according to 



- N 

IT a n f(A + /3 n ) 



7 



~T N \{[a n {A + l3 nl 1 - n ) 

n=l 



(16) 



where a n and (3 n follow according to Eqs. (|15fl from the 
parameters of Uqi and U%\ during the n'th passage and 
7 is given by AT = jTA, so that 



and 



A' 



F N {A) = J2 fnA n = f N [] (A + M 1 - 71 ) 



(17) 



(18) 



is an arbitrary polynomial of order TV in A determined 
by its roots — /3 n 7 1 ~". The product symbol Y[ nas been 
used to indicate that the index of the factors increases 
from right to left. 

Assume now that we want to implement an arbitrary 
expandable function F(A). In general, A can only be 
approximated by a polynomial F^(A) of order N in A 
with sufficiently large N. There are however two cases in 
which F(A) can be replaced exactly with some polyno- 
mial. 



In the first case, a set of N' + 1 eigenstates \<pi>) of A 
can be found. 



A\( Pl ,)=A' l ,\(p v ) = A l \( Pv ), 



1=0,. ..,N 
l'=0,...,N' 



(19) 



where N' (> N) is chosen such that a renumbering of the 
eigenvalues A[, gives a set {A , . . . , A^r} of exactly N + 1 
distinct A;. We now choose Eq. (|18fl to be the polynomial 



F N (A) 



N 

E 

1=0 



N 



f(aa n 



k=0 



A - A k 
At -At 



(20) 



of order N in A, so that we obtain Fn(A)\lpi>) ~ 
F{A)\(p v ). In Eq. (JTHJ) we then have 



where 



AT' 



fN 



■jk) _ 



(jit — ,3k) 
N' ) 



(21) 



JV' 

E 

l'=0 



(22) 



We see that a signal state g satisfying the condition 

~p(j'i.— jk) _ p(ii,— jk) ~ _ - 
(JJr N , — jt n , (J — Q 



(23) 



is transformed according to Eq. QJ, where Y ^T N F(A), 
cf. Eq. (|2*T|l . In this sense, a desired function F(A) can 
be implemented. In order to compensate the operator T, 
the transmittance T\ may be chosen sufficiently close to 
unity, i.e., \R\\ Cl. A precise compensation is possible 
by an additional implementation of the operator T~ N al- 
lowed in the photon-number truncated case as discussed 
below, in section IVI Al We then have Y ~ F(A). Note 
that the replacement of the function F(A) by a poly- 
nomial Fn(A) of order TV in A is unique. If there was 
another polynomial Gn(A) of order N in A for which 
also G N (A)\cpif) = F(A)\ifi>) with l' = 0,...,N',we could 
consider the polynomial Fn(x) — Gn(x) of order iV in 
x. It has however the (N + 1) distinct roots xi—Ai, so 
that Fn(x) — Gn(x) = 0. Of particular interest is the 
case in which in each of the modes j±, . . . ,jk, the number 
of annihilation operators occuring in A equals the num- 
ber of creation operators, so that the eigenstates \<pv) = 



are N' + 1 photon number states (they 
may be, e.g., the lowest states occuring in a Fock space 
expansion). 



B. Photon- number truncated signal states 

In the other case, there exists a mode ji in which the 
number of annihilation operators occuring in A exceeds 



5 



the number of creation operators and the signal state g 
is photon-number truncated in this mode, 

N 

Qp™ = P^g = g, P^ ] = £ In),-, A (n|. (24) 

We choose Eq. (|18f) to be the sum of the first N + 1 ele- 
ments of the power series expansion of F(A), so that the 
signal state g is again transformed according to Eq. ijTjl. 
whereY~f N F(A), cf. Eq. (ED- 



V. APPROACHING GENERAL OPERATOR 
FUNCTIONS 



VI. EXAMPLES OF APPLICATION 

We have seen that we may implement a given binomial 
in any moment of the creation and annihilation operators 
of a number of traveling optical modes. A repeated ap- 
plication then allows - in principle - the approximation 
of a desired function of the respective mode operators. 
In a concrete example, the general scheme can be sim- 
plified considerably. Let us illustrate this in some cases 
of special interest such as single- and two-mode opera- 
tors, since the application of single-photon doners may 
be avoided in this case. 



Assume that we want to implement an arbitrary func- 
tion F of the creation and annihilation operators of a 
number of signal modes, Y ~ F in Eq. QJ. We further 
assume that the expansion of In F can according to 



F 



JnF 



(25) 



be truncated after a (sufficiently large) number N of el- 
ements. Here, the A n are products of the respective cre- 
ation and annihilation operators that can be written in 
the form of Eq. I|15b|) . Consider now TV 2 successive ap- 
plications of Eq. (|14|l according to 

5tat = Y(A N 2 7 a N 2,l3 N ?) ■ ■ - Y(Ai,ai,f3i) 



- N 

TT \a n f{A n +f3„ 

AJ -n— 1 L 

- N 2 r 

J- A n— 1 . 



(26) 



where j n is given by A n T — j n TA n , cf. Eq. (|T7jl . Choos- 
ing A n as well as /3 n according to 



An 



,-1 ..n-1 



n = Nc^Y„ 
where [n] = 1 + [(n — 1) mod N], we see that 



N 



n=l 



(27a) 
(27b) 



(28) 



where 



X 



' N 



/V 



(iV>l) 



(29) 

and therefore X ivF for TV 3> 1. Again, we may compen- 
sate the operator T by choosing T\ sufficiently close to 
unity or subsequently implementing the operator T~ N , 
cf. section IVl Al so that we obtain Y to t ~F. 



Functions of the single-mode photon number 
operator 



Until now we have placed emphasis on the realization 
of general operators. In a given case however, the respec- 
tive scheme may be simplified considerably. For example, 
let us assume that A — ho. This example is of particu- 
lar importance since it allows the compensation of the 
operators T in Eq. and Eq. I|28[) for photon-number 
truncated states as we will see below. In Eq. (|14ll . we 
then have k — 2, j\ = = and s± = 0, S2 = 1. Instead 
of applying the setup Fig. [3 in the signal mode first in 
case (0) and then in case (1), it is sufficient to prepare its 
input ports 1 and 2 in a state \^)i2 as defined in Eq. J7J. 
The resulting complete setup is shown in Fig. Eta). If 
the photodetectors Di and D 2 detect and 1 photons, 
respectively, the signal state g is transformed in mode 
according to Eq. (JIJ, where 



Y(A,a,P) 



i(0| 2 (l|f/2iC/oi|*)i2 
alf°(A-f3), 



(30) 



with 



.4 



n , 



\/2Tx 
2 



Rx 



1 



T*R 2 



(31a) 
(31b) 

(31c) 



For fixed C/01 , P can be arbitrarily chosen by varying the 
parameters of U%i. As a consequence, there is no need 
for cloning devices. In order to realize a transformation 
of a photon-number truncated signal state [cf. Eq. Ij24(l 
with ji — 0] according to Eq. (JJJ with a desired function 
Y ~F(no), there have to be N successive applications of 
Eq. (|30|l with varying (3 n such that the overall operator 
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becomes 



Y tot = Y(A,a N ,f3 N )---Y(A,a lj f3 1 ) 

N 
n=l 



(32) 



where the symbol : : denotes normal ordering, and 
Eq. becomes 



JV 



v p(°) - 



Y[a n )N\(Tr N z-l)- J " z A °P^. (40) 



Inserting A — uq and with it Ai = I into Eq. (|2U|1 . and 
furthermore substituting T^~ Nl F(l) for F(l), we see that 
the (3 n in Eq. l].'{2[l must be chosen such that 

1 N fN\ N 
= Try E , ("2f )"'F(0 J] (fc-fio). (33) 



i=0 



fc=0 



Since then T^ fl ° F N (n )\n) = F(n )\n) holds for number 
states \n) with n = 0, . . . , N, we obtain 



JV 



itot-Pjv 



(0) 



[for 4 0) see Eq. ^2J] , where 



iV 



(_!)* " , N 
AT! ^ W 



(-Tf)-'F(Z), 



(34) 



(35) 



so that for photon-number truncated states the desired 
function F(no) is implemented. 

As an illustration, let us consider the special case of 
the exponential F(jiq) = z n °, which, if applied with ap- 
propriate z in modes jx,...,jh, allows the compensation 
of the operators T in Eq. i|ltjfl and Eq. I]28fl. Inserting 
F(l) — z l into Eq. (|33|) and applying 



fe=0 



x 



71 + 1 ( l)n( n 

x-r ' \l 1 Vn + l 



where x is a complex number, 1 = 0,1,..., and 
'x\ f 1 : A; = 



Hllto(z-j) : fc = l,2,. 



(36) 



(37) 



gives 



JV 



F w (n ) = £(Tf"z-l 



fe=0 
N 



fc=0 



A:! 



(38) 



Here, we have applied Eq. (|B2|) in the case of normal 
ordering, s = 1. Making use of Eq. (|B4(1 . we see that 



F N (n Q )P J 



fc=0 



. e (T 1 - JV Z -l)rl . 



p 



k\ ~ N 

(0) 



/V 



(39) 



B. Cross-Kerr nonlinearity 

Another example allowing a simplification of the gen- 
eral configuration is the implementation of a two-mode 
cross-Kerr interaction, Y ~ e lvnin ° . In Eq. ifTijl. we then 
have A = nin 0) fc = 4, ji = j 2 = 0, j 3 = j 4 = 1, and si 
= S3 = 0, S2 = S4 = 1. Instead of applying the setup 
Fig. n in the signal modes and 1 first in case (0) and 
then in case (1), a single application in each signal mode 
as shown in Fig. E^b) is sufficient. To the setup acting 
on mode 1 consisting of beam splitters Ui^(Tx,Ri) and 
^35(^2, R2) as well as photodetectors D 3 and D 5 , a sec- 
ond, identical setup acting on mode is added consisting 
of beam splitters Uo±(T\, R\) and Un(T2,Ri) as well as 
photodetectors D2 and D4. The modes are labeled as 
shown in Fig. |3Jb). The input modes are prepared in a 
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FIG. 3: Setups consisting of beam splitters £/ and photode- 
tectors D, conditionally implementing a desired binomial in 
no (a) and a desired binomial in (hi — 7) (no — 7) with some 
given 7 (b). 

state I 2345 as given by Eq. (|A9J| with modes relabeled 
accordingly. If the photodetectors D2 and D3 detect 1 
whereas D4 and D 5 detect photons, respectively, the 
signal state g is transformed in modes and 1 according 
to Eq. ||TJ, where 



Y(A,a,0) = 23(1|4 5 (0|[/24?704C/35^15|*)2345 

= a7f 1+ "°(i-/?), (41) 



7 



with 



A = (iii- 7)("o ~ 7), 



1 2 







7 



1 

Ti 



12 

#2 



(42a) 
(42b) 

(42c) 
(42d) 



splitter operator Ujk defined in the beginning by Eq. ((3J), 
which can equivalently be written in the form of |43|. 



U jk = (VTpe 



-P-R.*ala je V*Ka]a 



>(P*Ty 



(46) 



For fixed £, T\ and R\, (3 can be chosen as desired 
by varying T2 and R 2 . By N repeated applications of 
Eq. I|41|) with varying (3 n we obtain 



Its parameters are here written in calligraphic style to 
distinguish from those of the beam splitters in the de- 
vice. The implementation of exponentials of S^fij with 
k = j explained in section IVI Al represents the case of 
section IIV Al and those with k 7^ j represents the case of 
section FtV Bl Together with a detection of the vacuum 
state, the latter can be used to achieve a teleportation of 
a state g from mode j to mode k, 



Ftot = Y(A,a N ,/3 N )---Y(A,a 1 ,p 1 ) 



N 



where 



n /* t i 



lrpN(hi+ho) 



F, 



(43) 



\n=l 



We now tune T\ and R\ such that r ) = Lp~ 1 2kir with some 
arbitrary fc = 1, 2, • • • . The f3 n are chosen such that in 
Eq. 143|) the polynomial TV of order N in A becomes 



Qk — IkjQjljki 



n=0 



(47) 



(48) 



N 

1=0 



N 

n 

k=0 
(fe#i) 



A - A k 
Ai - A k 



(44) 



where Aq, . . . , Am are the (N + 1) distinct values of the 
expression (n\ — 7) (no — 7) with no,i = 0, . . . , M [cf. the 
remarks following Eq. 119|) ]. On signal states g whose 
photon number is limited to some given value M in mode 



A comparison with Eq. i|46|) shows that exp(aJ,a J ) [or 
exp(a|a fc )] is realized for T = 0. Physically, the beam 
splitter then describes a mirror. Since Eq. (|46|l becomes 
singular in this case, we have to construct it by a suc- 
cessive implementation of two operators Eq. 146|) whose 
combined action results according to 

Ujk{T 2 ,n 2 ,'p 2 )u jk (T l ,K 1 ,Vi) 
= U ih (r 2 T x -KJl\,TJl x +ll^,V 2 V x ) (49) 



and 1, i.e., g satisfies Eq. (gljl, where P^ l> is replaced again in a U(2) beam splitter operator. Choosing T 3 = 

IZj = 1/V2 and Vj = 1, we get an operator Ujk(0, 1, 1). 



with Pjf? ® , Eq. (|43|) then acts like an operator 



' iV 

n 



f — 1 N (fix +n ) iipA 1 n 
JN 1 \ e 



(45) 



In order to implement the unitary operator e 11 ^™ 1 ™ , Ti 
may be chosen sufficiently close to unity. Again, there 
is the alternative of a precise compensation of the extra 
exponential in Eq. I|45|) by an additional implementation 
of the operator j , ~ Ar (" 1+n °) as discussed in the previous 
section IVI Al In this way, we have proposed a beam 
splitter arrangement relying on zero- and single-photon 
preparations and detections that acts on photon-number 
truncated states like a two-mode Kerr nonlinearity. 



C. U(iV) 

Other operators whose action on photon-number trun- 
cated signal states can be realized exactly by a polyno- 
mial are those performing a U(iV) transformation. Since 
they can be factorized into U(2) couplers Ujk of the re- 
spective modes, let us limit attention to the U(2) beam 



If implemented by our scheme with separated modes j 
and k, such a 'quantum mirror' realizes a bidirectional 
teleportation between these modes. 

More generally, we may consider a transformation in 
N modes j = 1, . . . , N of a signal according to 0) by an 
operator Y ~U whose action is defined by 



(50) 



where p(l), ■ ■ ■ ,p(N) is a permutation of 1, . . . ,N. The 
set of such permutations of N mode indices forms the dis- 
crete symmetric subgroup S(N) of U(iV) physically rep- 
resenting iV-mode mirrors (with the identical operation 
included). Note that in general we may add additional 

phase shifts according to U''=e 1 ^'^= 1 VjUl U, so that S(N) 
is replaced with U(1) JV " <£> S(N). If U is implemented by 
the setup Fig. |2 with spatially separated signal modes, 
then we again have the situation of a teleportation. One 
possibility to construct U is to introduce auxiliary modes 
— 1, . . . , — N, and to apply l|4*5|) repeatedly according to 



U 



1 p{N),-N 



■L 



p(i), 



-N,N ■ 



(51) 
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D. Preparation of multi-photon entangled states 

We conclude with an example of state preparation 
which may be understood as a state transformation with 
a given input state. Consider the following situation. By 
feeding the overall setup with vacuum, g = |0)(0|, the 
scheme may be used to generate a state difficult to pre- 
pare by other means. An example are fc-mode states 



\z)i,...,k 



l-\z\ 



N 



1 - \z\ 2 ( N +V ^~ 

1 1 71=0 



(52) 



which can be generated by applying a polynomial 



Y 



N 

E 1 

n=0 



(zA) n 



(53) 



in A= a' k ■ ■ ■ a[ to the vacuum. The states Eq. (|52(l are of 
interest from the theoretical point of view. Consider the 
case \z \ < 1 and N^oo. We see that for k = 1 we obtain 
a coherent phase state (in particular, for \z\ = 1 a London 
phase state). For k — 2 we obtain a two-mode squeezed 
state (in particular, for \z\ = 1 an EPR-like state). In 
the limit k — > oo, the states Eq. 13211 may be used as a 
representation of an EPR-like state between a number n 
of quantum fields (instead of two modes of a field) . To 
see this, the fields are indexed according to = 1 + [(j 
— 1) modn], i.e., mode j is ascribed to field . 



VII. 



CONCLUSION 



In summary it can be said that, starting from a sin- 
gle photon, cf. Eq. Q, we can with a given probability 
approach arbitrary operators of spatially separated trav- 
eling optical modes. The problem is that the success 
probability, which depends on the desired transforma- 
tion and the signal state in general, is expected to de- 
cline exponentially with an increasing number of steps 
N, so that the practical applicability of the scheme is 
limited to small N, sufficient to engineer, e.g., states in 
the vicinity of the vacuum. On the other hand, the ap- 
plication of giant nonlinearities to fields containing only 
few photons represents just the case discussed in poten- 
tial quantum information processing devices. To give an 
estimation, consider a passively mode-locked laser whose 
resonator has an optical round trip length of I = 3 cm, so 
that the repetition frequency of the emitted pulse train 
is v = l~ x c = 10 GHz. If this pulse train is used to 
prepare the entangled state needed to run the scheme, 
then its repetition frequency determines that of the over- 
all device. An assumed total success probability of p = 
10~ 10 would then reduce the mean repetition frequency 
of the properly transformed output pulses to v' = 1 Hz, 
which may still be useful for basic research experiments 



but unacceptable for technical applications. Apart from 
this, there may be the situation where only a single un- 
known signal pulse is available that must not be spoiled 
by a wrong detection since no copies can be produced of 
it. This constitutes the main drawback of all conditional 
measurement schemes. 

It may be worth mentioning that the whole scheme 
may also run 'time-reversed', i.e., all pulses are sent in 
opposite directions through the device, while the loca- 
tions of state preparations and (assumed) detections are 
interchanged, so that the resulting 'adjoint' scheme im- 
plements the adjoint operator YK 
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APPENDIX A: SINGLE-PHOTON CLONING 

In what follows, we give a number of possibilities to 
implement a single-photon doner Eq. JSJ. 

1. Cross-Kerr nonlinearity 

To implement single-photon cloning, we may apply a 
Mach-Zehnder-interferometer equipped with a cross-Kerr 
nonlinearity as can be seen as follows. We prepare input 
mode 2 of a balanced beam splitter U21 (T = — R = 1/ v2) 
in a single-photon state |1) and the other input mode in 
the vacuum state |0). After output mode 2 has passed a 
cross-Kerr coupler 



U 



(Al) 



it is remixed with output mode 1 at a second beam split- 
ter U^i- The output mode 2 of this second beam splitter 
is then again mixed with an auxiliary mode 3 prepared in 
a coherent state |£) with £ = T _1 i? using a beam splitter 
U32 (T, R) . If eventually and 1 photons are detected in 
modes 2 and 3, respectively, the action of the complete 
setup on the state of the input mode can be described 
by an operator 



Y 



2(0| 3 (l|[/ 32 f7 2 t 1 f/o2?72l|0)l|l) 2 |e)3 
_ 1 . 1 _|_ f_1 "itio+s 

3 iarsfl ^E i±i ^ — i->i. 

s=0 



so that 



YR 



(0) 



_i are R 



Vp Qi 



[for Pp see Eq. iJUJl]. The success probability 
P= f(l + l£|- 2 )e lf|2n ~ 



(A2) 



(A3) 



(A4) 
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attains for |£| 2 = (a/5 - l)/2 w0.62 its maximum value of 

Pmax~0.21. 



Eq. (|A10|) we have relabeled output mode 2 as mode 0. 
The success probability 



2. Three- wave-mixer 

A disadvantage of the above proposal is its dependency 
on a given measurement result. Unitary photon cloning 
can be achieved using a three-wave-mixer 



#012 (<p) 



= pi¥>(a a i a 2+ a o a i a 2) 



(A5) 



Preparing its input mode in a photon number state \s) 
with s = 0, 1 and its input modes 1 and 2 in the vacuum 
state |0), the output state becomes 

#oi2(y)|s)o|0)i2 = cos(s^)|s) |0)i2 + isin(s<p)|0)o|s)i 2 , 

so that for a phase of ip — ir/2 and an additional appli- 
cation of a preceding phase shifter e -1 ^™ we obtain an 
output state Y\s)o = |s)i2, where 



Y = o(0|#oi2(7r/2)e-^ fto |0)i 2 . 
After relabeling output mode 2 as we get 



YR 



(0) 



(A7) 



(A8) 



cf. Eq. JSJ) and Eq. I|24|l . Formally, we could consider a 
success probability for which p = o(s\Y'Y\s)o = 1 holds 
as expected for an unconditional unitary operation. 



3. Linear optics 

A disadvantage of the previous two proposals is their 
dependency on large nonlinear coefficients which are 
hardly achievable in practice. Therefore, we now give a 
possibility of implementing Eq. JSJ) solely based on beam 
splitters and zero- and single-photon detections. Assume 
that a state 



I*) 



1234 



|l)l2|0)34+g|0)l2|l) 3 4 



(A9) 



with some arbitrary £ ^ is prepared. Mode 3 is mixed 
with an auxiliary mode 5 prepared in a coherent state 
|£) using a beam splitter U 53 (T = R= l/v2) and mode 4 
passes a beam splitter U ±(T = R = 1/a/2)- If eventually 
photons are detected in modes 3 and 4 but 1 photon 
in each of the modes and 5, respectively, the device 
realizes an operator 



Y = 34<0| 05<l|#04#53|*)l234|0l 



3 i arg£ 



VpQi 



(A10) 



that acts on states \s) in which input port of the beam 
splitter £/q4 is prepared. Note that in the second line of 



4(i + ir> 



(All) 



becomes maximum for |£| 2 = (a/5- l)/2«0.62 for which 
p«0.05. 

One way to prepare the state Eq. (|A9|) from single- 
photon states is shown in Fig. 0] If transmittance and 
reflectance of the beam splitter U Qi are denoted by T and 
R, respectively, then those of beam splitter C/35 are given 
by T* and R* . The parameters of all other beam splitters 
are given by T = i?=l/\/2. The setup represents an array 



|i) 1 



11} < 5 



|1) ,0 



11} X 2 




FIG. 4: Setup consisting of beam splitters U and photode- 
tectors D, conditionally preparing a state given by Eq. fAlil. 

of these beam splitters which is fed with a number of 
single-photon states |1) and vacuum states |0) as shown 
in the figure. If each of the two photodetectors detects 1 
photon, the output state becomes 



^05(l|^50^34C/o4t/35^12[/ 02 c/i5|l)oi25|0)34 = |*)l234, 

Vp 

( Al2 ) 

where ip = tt — 2 arg T is an unrelevant phase and 



6 = -R* 

The success probability reads 



P 



i-iei 



(1 + iei 2 ). 



(A13) 



(A14) 



We see that p drops from its maximum value of 0.25 as 
attained for £ = to if |£| — > 1. For the above value of 
|£| 2 = 0.62 we obtain p« 0.02. 



10 



APPENDIX B: s-ORDERING RELATIONS FOR 
THE PHOTON NUMBER OPERATOR 

In what follows, we compile some relations regarding 
the s-ordering of the photon number operator n = a) a 
as used in section IVI Al We start with its fcth power. 
Introducing s-ordering according to |44j . 



minima) , . . . /+_,\ k 



k=0 



and applying [45j, p. 626, no. 30, we get 



(Bl) 



' 1 — s\ k 1 fk\ fii 



1=0 



I J \ l 



(B2) 



The inverse relation is obtained from Eq. I|B1|) with s = 
1 and t renamed as s, inserting Eq. (|B2fl with s = 1, and 
applying ^5|, p. 624, no. 25, as well as p. 614, no. 31, 
which yields 



To consider the case where h is the exponent, we apply 
HU, p. 709, no. 3, as well as p. 612, no. 1, to Eq. (p2jl . 
which gives 



(i 



The inverse relation 



(1 + 



n+l 



1 I (1 ~ S) 2 ( "~~ 



(B4) 



(B5) 



m,j,l— 



s— 1 \ u 



E ( " )(-l) m mM'^VL ' 2 " 



mj-0 



1-s 



then follows directly from Eq. (|B4|) in accordance with 
.(B3) 0. 
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